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Abstract 

A beam of diatomic molecules scattered off a standing wave laser mode 
splits according to the rovibrational quantum state of the molecules. Our 
numerical calculation shows that single state resolution can be achieved 
by properly tuned, monochromatic light. The proposed scheme allows for 
selecting non-vibrating and non-rotating molecules from a thermal beam, 
implementing a laser Maxwell's demon to prepare a rovibrationaUy cold 
molecular ensemble. 
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Manipulating the motion of massive particles by using the mechanical ef- 
fect of light has been a rapidly developing field since the discovery of lasers. 
Understanding the fine details of the light-matter interaction on a microscopic 
level contributed much to the advent of impressive applications, including laser 
cooling and trapping of atoms and atom optics Q]. The adaptation of these 
techniques for molecules would undoubtedly open new perspectives. Recently, 
pioneering experiments have been reported in this direction, for trapping ^ 
and decelerating [Q molecules or for interferometry [|j , and also for fundamental 
tests of quantum mechanics |^ with molecules. It is, however, nontrivial to ex- 
tend the optical methods to molecules because of the complexity of their energy 
level structure. Optical cooling of molecules would require new approaches, 
for example, such as coupling the particle to a far-detuned light field inside an 
optical cavity (s). 

The essential difference stems from the additional mechanical degrees of 
freedom. The vibrational and rotational modes lead to densely spaced levels. 
Typical excitations are orders of magnitude smaller than the electronic transi- 
tions. A usual thermal beam consists then of molecules in various rovibrational 
quantum states. Cooling translation, vibration, and rotation simultaneously is 
a difficult problem; a complex scheme with multiple laser beams has been pro- 
posed in 1^. However, there is another option to obtain molecules with lower 
energy: to sort the thermal ensemble by means of state-sensitive interactions. 
A Stern-Gerlach-type experiment with inhomogeneous microwave field has been 
first performed to separate the rotational states constituting the resonant tran- 
sition within a given vibrational band jl^ . Rotational state dispersion has been 
shown to be present in various molecular optics settings using far-detuned laser 
fields . In this paper we show that a Maxwell demon like operation can be 
achieved in a different regime of the molecule-light interaction. That is, a single 
rovibrational state, possibly the ground state, can be filtered from a thermal 
beam. Owing to the high degree of control facilities, the proposed scheme can 
be applied to a quite general class of molecular states. 

Deflection of a CS2 and I2 beam using the nonresonant dipole force has 
recently been demonstrated fl^ . The laser was well detuned from a given 
electronic transition and all the molecules experienced very nearly the same 
nonresonant dipole force regardless of what their internal quantum state was. 
In order to accommodate the nonresonance condition for all the transitions, 
detuning on the scale of the electronic transition frequency was applied. More 
sophisticated control of the mechanical effect of light could be obtained with 
a laser whose bandwidth is comparable with the linewidth of the rovibrational 
transitions in the spectrum. With detunings on the order of the rotational 
level-spacing, less laser power is needed to have the same mechanical effect. 

The reaction of atoms and molecules to radiation is usually described by 
the frequency-dependent polarizability tensor. This is a good approximation 
for a transparent medium, when the radiation frequency is far away from any 
resonances in the matter. The polarizability can then be calculated from the 
quantum mechanical mean value of the dipole operator, usually by employ- 
ing perturbation theory. This approach is sufficient to interpret the deflection 
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Figure 1: Scheme of deflecting molecules quantum state selectively. A beam of 
diatomic molecules with velocity vq is scattered off a monochromatic, standing 
wave laser mode. The detector D measures the deflected beam intensity along 
the z axis. 

experiment and many other schemes for manipulating molecular motion 
[ p^ [l^ . However, the simple concept of a state-independent polarizability does 
not hold when the monochromatic laser frequency is close to some individual 
rovibronic resonances. In this case one must return to a more fundamental 
quantum mechanical treatment relying on the dipole operator itself. Otherwise 
the standard theoretical framework for treating molecular optics ||T^ , that is the 
adiabatic separation of the rovibrational modes and the center-of-mass motion 
can be used. 

We consider a beam of diatomic molecules (see Fig. 1) moving with ve- 
locity Vq in the x direction and crossing at a right angle a CW standing- 
wave laser mode. The interaction energy, in the dipole approximation, reads 
V = — dE. The electric field is linearly polarized, and its spatial variation is 
given by E = Eo cxp [— (x^ -I- y'^)/2w'^] cos kz cos cut. The molecules experience 
a force proportional to the gradient of the inhomogeneous interaction energy. 
Due to this light force, the molecules accumulate transverse momentum, Mv^. 
However, the displacement in this direction is negligible during the interaction 
(Raman-Nath regime). The laser-induced potential is assumed to be much 
smaller than the kinetic energy; hence the longitudinal velocity = vq remains 
constant during the whole passage. 

The spectrum of the diatomic molecule in the Hund's case a is defined |^ 
by the Hamiltonian 

Ho=Y^ E{n, V, J, M, n) \nvJMM:) {niyJAIAY.] , (1) 

where n labels the electronic states, v is the vibrational, J and M are the usual 
rotational quantum numbers, A and S are the projections of the total electronic 
angular momentum and spin on the intermolecular axis, respectively, and f2 is 
their sum. Nuclear spin is ignored; however, we will comment on this point later 
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on. Within the Born-Oppenheimer approximation, the electronic and nuclear 
wave functions can be separated as follows 



\nvJMhYs) ^ <S>nA^(Te,r)RnAr)ujMn{0,(f>) , (2) 

where the variable Te incorporates all the electron coordinates, r is the internu- 
clear distance, and 6 and (p define the orientation of the molecule. The electronic 
spin part of the quantum state is omitted. The energies for low v and J are 
usually approximately given by the expression 

E{ n , J, M, n) = 

E,i{n) + huJe{y + 1/2) - huJeXe{f + 1/2)^ 

+B^[j{j + 1) - n^] - D^[j{j + 1) - n^f , (3) 

where the first term is the electronic energy, the next two terms describe the 
anharmonic vibration, and the last two correspond to the energy of a symmetric 
top (the standard notation is used). 

Let us suppose that only the fundamental / and one excited e electronic 
states are effectively coupled by the electromagnetic mode. The interaction 
Hamiltonian can be restricted to 

•E E Ar::fM'\f^JM){ev'J'M'\+h.c. (4) 

uJM u'J'M' 

The electric field Eq is linearly polarized in an arbitrary direction perpendicular 
to the mode axis z. It is enough to consider the dipole moment in the direction 
of the polarization. We resort to the Franck-Condon approximation to calculate 
the matrix elements of the dipole operator, i.e., the induced electric dipole mo- 
ment is supposed to be independent of the intcrnuclcar distance. This separation 
is not fundamental, but only convenient for purposes of discussion. The matrix 
element can then be decomposed into three factors, rfj^/jv/^ = d'jJi'^ ^jMn^ ■ 
The electric dipole is given by the integral 



d}^ y" (Te , r)d$ / (Te , r)dT, , (5) 

averaged over the nuclear distances r. The next term, 

r: = [ R:,,Rfydr , (6) 



is the vibrational overlap factor (the square root of the Franck-Condon factor) 
that we will calculate using Morse-potential eigenfunctions. Finally, a rotation- 
dependent factor follows which can be analytically evaluated for the symmetric 
top eigenfunctions. Two cases can be distinguished: 



-j'M'n' 



COS 

u*j'M'n'UjMn { gjj^p \sm9d0d(l), (7) 
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where cos 9 is used for transitions with Afl = 0, and sin6' is used for transitions 
with AO = ±1. 

The Hamiltonian, in the rotating wave approximation, can be transformed 
to the form 

^- -E E 

lyJM u'J'M' 

KAf {\f^JM) {fvJM\ + \ev'J'M') {ev'J'M'\) 

+ ileiy'J'M') {eu'J'M'\ - \fiyJM) {fvJM\) 

+ng:fij'''f{Il) {\fvJM) {ev'J'M'\ + h. c.) . (8) 

The first term defines the energy zero level for each rovibronic transition \ fvJM) < 
\ev' J' M') . The second term is the corresponding detuning from the laser fre- 
quency, Mlf^' = E{ev'J'M') - E{fvJM) - fiu. The coupling strength is 
given by g^jM^ = ^Y^jm^ Eo/?i. The mode function /(R) collects the spatially 
varying factors. The coupling constant can be expressed in terms of measurable 
quantities as 

where Ty-^/ is the natural linewidth, and S{J,J') is the Honl-London factor. 
The laser parameters are the intensity / in units of photons/second and the 
effective beam area A. 

The fine tuning of the laser is the crucial point of our scheme. The laser in- 
duces nonresonant interactions if the condition {5^jm ^ )^ ^ idvjM ^ )^ holds for 
each transition. The laser is detuned from all the transitions, > Y'^jm ^ , 

and the smallest detunings are on the order of the rotational level spacing B^. 
The nonresonance condition on this scale limits the coupling g^jjj^^ <C B^, and 
thus the laser intensity. 

The molecules are initially in the electronic ground state, and the vibrational 
and rotational states are populated according to a thermal distribution. The 
Hamiltonian (H) describes the laser-induced dynamics of the molecule except for 
the translational motion. Interaction with the vacuum modes yielding sponta- 
neous emission was omitted, since the upper electronic states are not populated 
by the nonresonant excitation. The translational motion is treated classically in 
our model. The center-of-mass (CM) coordinate appears as a parameter in H^f^ 
via the mode function /(R). Motion in the y direction is irrelevant with respect 
to the deflection scheme. We insert v^t in the coordinate x, thus making the 
problem explicitly time-dependent. For slow enough motion, the system evolves 
adiabatically. In the nonresonant limit the eigenstates of H^tt, which can be 
viewed as dressed rovibronic states, are close to those of the unperturbed Hq. 
Transitions can be neglected between these eigenstates. The eigenenergies are 
shifted (AC Stark effect) by an amount that can be well approximated by taking 
only the most significant transition into account for every rovibronic state. In 
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fact, this approximation is equivalent to the decomposition of H^tt into 2 by 
2 blocks. Corrections can be generated systematically by extending the block 
size with the inclusion of weaker and weaker couplings. The consistency of this 
method is based on the hierarchy among the couplings that can be set up ac- 
cording to the g^jM ^ /Kjai^ ratios. This is not a perturbational expansion: 
whenever an additional coupling term in H^tf is considered, its contribution is 
completely included. In a first approximation, the energy associated with the 
state \ fvJM) reads 



^jfir ± h^{g:'£rYf{w + {K'uryi'^ , (lo) 

where E'^/m^ , g^jht^ ■> ^-nd are the parameters of the transition to the 

strongest coupled upper rovibronic state \ev' J' M') . The sign is determined by 
the sign of the detuning. During the passage the molecules acquire transverse 
velocity due to the nonresonant dipole force, given by the gradient of the en- 
ergy Eq. (pX|). We predict the occurrence of different trajectories labelled by 
the quantum numbers of the relevant rovibrational states. The quantum state 
dependent deflection angles read 

a{iyJM) j^, , (11) 

"x "uJM 

where I = \/A is the effective interaction length, = hk/M is the photon 
recoil velocity. 

In the following, we analyze the complex deflection pattern on the example 
of the Na2. The ground electronic state X^T,^ is coupled to the state ^^S^ by 
a dipole-allowed optical transition with E^i = 17541 cm^^ |l^. Both states are 
singlet; hence the electron spin plays no role (E = E' = 0). Since A = A' = 0, the 
selection rules J' — J ± 1, A/' = M apply (the rotational Q branch is missing) 
and the Honl-London factor is S{J,J') — J + 1. According to theoretical 
calculations |l^, confirmed by experiments |2^, the electric dipole is |djp = 14 
in atomic units. 

Suppose the laser is tuned close to resonance with the rovibronic transi- 
tion /, z/ = 0, J = 0, M = ^ e, I/' = 6, J' = 1, M' = (with energy 

+ 666.136 cm~^, and natural linewidth T = 3.4 x 10^^ cm^^). The ^ 6 vi- 
brational transition is associated with one of the largest Franck-Condon factors 
(|i?Qp = 0.113). The neighboring transitions are represented in Fig. 2a. For a 
temperature of 1000 K, the distribution can be truncated at 10 vibrational and 
100 rotational levels. As the coupling strength is much less than B^, resonances 
that are as far as the rotational line separation {B^ = 0.155 cm"'^) can be safely 
neglected. Thus a molecule has a chance to be deflected if its rovibrational state 
is involved in the transitions shown in Fig. 2a. The adjacent lines are separated 
by at least 0.05 cm~^; hence g « 1.5 x 10~^ cm~^ can be chosen (F <C (? <C B^). 
Focusing the standing wave to / ~ 50/im, the necessary field intensity is / « 0.3 
mW. The Raman-Nath regime implies a purely geometrical condition on the 
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Figure 2: a) The spectrum of Na2 around the transition /, 1^=0, J=0, M=0 
<-> e, i/'=6, J'=l, M'=0. E^i is subtracted from the frequencies, b) Deflection 
angles as a function of the laser frequency. Dashed vertical lines correspond to 



the transitions. Outside the grey regions 5^ 



v'j'M' 
uJM 



> 10 



v' J'M' 
vJM 



thus the non- 



resonance condition is fulfilled. At the specific laser tuning loq several deflection 
angles occur (points). The arrow indicates the deflection of molecules in the 
vibrational-rotational ground state. 
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deflection angle a that is limited by a < A/Z ~ 10~^ radian. The velocity is 
set at = 500 m/s, yielding an interaction time of 100 ns. The probability of 
a spontaneous emission event is less than 10 %; hence neglecting the scatter- 
ing force is justified. Deflection angles according to the above parameters are 
shown in Fig. 2b as a function of the laser frequency. Each curve belongs to a 
single rotational-vibrational state. The resonant region around each transition 
is darkened. Outside these regions the detuning is one order of magnitude larger 
than the corresponding coupling constant, thus the nonresonance condition is 
fulfilled. For the laser frequency cjq at the solid vertical line, the detuning will be 
5 K, 0.02 cm~^. The points represent the observable deflection angles. The ar- 
row indicates the angle corresponding to the rotational-vibrational ground state 
(approximately 100 /irad; for comparison, a single photon recoil would cause 35 
^rad). According to Fig. 2b, several deflected partial beams are present simul- 
taneously, which could be resolved depending on the beam divergence and on 
the detection technique. We note that a deflected beam does not correspond 
to a single quantum state, as we did not include in our model all degrees of 
freedom of a molecule, e.g. the nuclear spin was neglected. Thus the deflected 
beam contains all possible nuclear spins, altogether 6 states. 

The center-of-mass motion has been described in terms of classical trajec- 
tories in our model, which offered an intuitive picture with focus on the more 
essential internal part of the dynamics. We avoided, however, some details that 
would degrade the predicted signal in an experimental realization. There are 
several mechanisms giving rise to a broadening of the deflected partial beams. 
First of all, we assumed a well-defined longitudinal velocity Vx ■ It follows from 
( pT| ) that the double the relative variance in Vx equals the relative variance 
of the deflection angle. Still we can safely tolerate a few percent uncertainty 
which can easily be achieved, even without velocity selection, by using super- 
sonic beams. Next, the molecule has to cross the field in a region where the 
interaction strength, i.e. /(R), is approximately uniform. The defiection angle 
being proportional to sin2fcz, the crucial part in this respect is the direction 
z. The variation 2fcAz — ±0.45 around 2kz = tt/2 amounts to 10 % varia- 
tion in the final deflection angle. However, the confinement of the beam to 
a size smaller than the wavelength causes diffraction. The transverse velocity 
that was supposed to be fz = broadens to a distribution with a width of 
?i/(2MAz) = w„<,/(2fcAz) « v,.^^. Deviations in the direction y from the center 
line y — Q do not induce such drastic diffraction since the gaussian beam waist 
is several orders of magnitude larger than the wavelength. Therefore « 
remains a good approximation. Finally, the spontaneous photon emission is 
accompanied by a recoil that introduces a randomness in the deflection angle. 
When the spontaneous emission probability is kept below 10 %, the occurence of 
several photon recoil kicks can practically be excluded. Similarly to the diffrac- 
tion problem in the z direction, the single photon recoil raises the question of 
resolution. That is, the deflection broadening, at most v,^^/vx, must be domi- 
nated by the deflection angle given in (^ij). Note that these two quantities were 
35 /irad and 100 /irad in our numerical example. 

The effect of spontaneous emission is, however, a bit more complex than 
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just broadening the velocity distribution. Following a spontaneous jump, the 
molecule exits the two states of the quasiresonant transition and the process of 
accumulating transverse momentum is interrupted. This gives rise to a weak 
background between the zero and maximum deflection angles. The suppres- 
sion of this background to the possible minimum is the guiding principle when 
chosing the appropriate detunings. In fact, this imposes the ultimate condition 
on the generalization of the proposed scheme to other molecules. The basic 
ingredient of our scheme is the "isolated line" in the spectrum. That is the 
spectral density must leave enough space to tune the laser quasiresonant with 
a single transition without introducing photon scattering. Roughly, an order of 
magnitude between the detuning and the linewidth V^fjJ^ , and also 

an order of magnitude between and the distance to the adjacent transi- 

tions in the spectrum is necessary. In Na2, being a light diatomic molecule, we 
had three orders of magnitude difference between T and the relevant transition 
line separations. This suggests that the non-resonance condition can also be 
fulfilled for heavier diatomic or polyatomic molecules as long as the rotational 
constant B^, is not less than one tenth of that of the sodium we considered here. 
Decisive conclusions on the applicability cannot be drawn without the detailed 
knowledge of the spectrum. Nevertheless, it may always be possible to find 
accidentally well isolated lines associated with the ground rovibrational state of 
the molecule. 

In conclusion, we have shown that it is possible to deflect diatomic molecules 
quantum state selectively. The very high sensitivity originates from the monochro- 
matic laser field that can be finely tuned within the rovibrational band. The 
laser-induced coupling of the vibrational and rotational modes to the center-of- 
mass translation results in quantum correlations between the different motional 
degrees of freedom. No analogue of this type of correlations exists in atom 
optics. There only the quantized light mode (quantum prism pl[|) or the in- 
ternal electronic state (optical Stern-Gerlach effect |2^) can be entangled to 
the translational motion. By contrast to atoms, however, in a usual beam of 
molecules various rovibrational states are thermally populated; hence the pos- 
sibility of selecting the rovibrationally cold molecules has important practical 
consequences. This arrangement may be integrated in the preparation stage 
of experiments performing coherent manipulation with molecular beams. The 
process still cannot be called cooling, since thermal equilibrium is not main- 
tained throughout. In fact, it is a laser Maxwell's demon which is not in 
contradiction to the second law of thermodynamics if the momentum exchange 
with the laser is correctly taken into account The Maxwell's demon for 

molecules, beyond its unusual statistical properties, would be of great practical 
value for preparing a sample of non- vibrating and non-rotating molecules. 

We are grateful for useful discussions with Jorg Schmiedmayer and Viktor 
Szalay. This work was supported by the National Scientific Fund of Hungary 
(OTKA) under contracts No. T023777, F032341, and F032346. T. K. acknowl- 
edges the support of the Hungarian Academy of Sciences (Bolyai Janos Kutatasi 
Osztondfj). 



9 



References 

[1] S. Chu, Rev. Mod. Phys. 70, 685 (1998); C. Cohen-Tannoudji, ibid. 707 
(1998); W. D. Phillips, ibid. 721 (1998). 

[2] T. Takekoshi, B. M. Patterson, and R. J. Knize, Phys. Rev. Lett. 81, 5105 

(1998). 

[3] J. D. Weinstein, R. deCarvalho, T. Guillet, B. Friedrich, and J. M. Doyle, 
Nature 395, 148 (1998). 

[4] H. L. Bethlem, G. Berden, and G. Meijer, Phys. Rev. Lett. 83, 1558 (1999). 

[5] M. S. Chapman, C. R. Ekstrom, T. D. Hammond, R. A. Rubinstein, 
J. Schmiedmayer, S. Wehinger, and D. E. Pritchard, Phys. Rev. Lett. 74, 
4783 (1995). 

[6] M. Arndt, O. Nairz, J. Vos-Andreae, C. Keller, G. van der Zouw, and 

A. Zeilingcr, Nature 401, 680 (1999). 

[7] For a review see J. T. Bahns, P. L. Gould, and W. C. Stwalley, Adv. At. 
Mol. Opt. Phys. 42, 171 (2000). 

[8] V. Vuletic and S. Chu, Phys. Rev. Lett. 84, 3787 (2000). 

[9] J. T. Bahns, W. C. Stwalley, and P. L. Gould, J. Chem. Phys. 104, 9689 
(1996). 

[10] R. M. Hill and T. F. Gallagher, Phys. Rev. A 12, 451 (1975). 
[11] T. Seideman, J. Chem. Phys 107, 10420 (1997). 

[12] H. Stapelfeldt, H. Sakai, E. Constant, and P. B. Corkum, Phys. Rev. Lett. 
79, 2787 (1997); H. Sakai et al, Phys. Rev. A 57, 2794 (1998). 

[13] J. Karezmarek, J. Wright, P. Corkum, and M. Ivanov, Phys. Rev. Lett. 82, 

3420 (1999). 

[14] B. Friedrich and D. Herschbach, Phys. Rev. Lett. 74, 4623 (1995). 

[15] T. Seideman, J. Chem. Phys. 106, 2881 (1997); T. Seideman, ibid. Ill, 
4397 (1999). 

[16] I. Kovacs, Rotational structure in the spectra of diatomic molecules, 1969, 
Akademiai Kiado (Budapest) & Adam Hilger Ltd (London). 

[17] G. Herzberg, Molecular spectra and molecular structure. I. Spectra of di- 
atomic molecules, 1950, D. van Nostrand Co. (Toronto). 

[18] B. M. Smirnov and A. S. Yatsenko, Physics Uspekhi 39, 211 (1996). 

[19] W. J. Stevens, M. M. Hessel, P. J. Bertoncini, and A. C. Wahl, J. Chem. 
Phys. 66, 1477 (1977). 



10 



[20] T. W. Ducas, M. G. Littman, M. L. Zimmerman, and D. Kleppner, J. 
Chcm. Phys. 65, 842 (1976). 

[21] P. Domokos, P. Adam, J. Janszky, and A. Zeilinger, Phys. Rev. Lett. 77, 
1663 (1996). 

[22] T. Sleator, T. Pfau, V. Balykin, O. Carnal, and J. Mlynck, Phys. Rev. 
Lett. 68, 1996 (1992). 

[23] V. S. Letokhov, Contemp. Phys. 36, 235 (1995). 

[24] G. J. Milburn, Austrahan J. Phys 51, 1 (1998). 



11 



